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We show that opacity of a clean multilayer graphene flake depends on the helicity of the circu- 
lar polarized electromagnetic radiation. The effect can be understood in terms of the pseudospin 
selection rules for the interband optical transitions in the presence of exchange electron-electron 
interactions which alter the pseudospin texture in momentum space. The interactions described 
within a semi-analytical Hartree-Fock approach lead to the formation of the topologically different 
broken-symmetry states characterized by Chern numbers and zero-field anomalous Hall conductiv- 
ities. 



I. INTRODUCTION 

Multilayer graphene is a link between one atom thick 
carbon layers^ with the peculiar Dirac-like effective 
Hamiltonian for carriers^ and graphite^ which can be 
seen as millions of graphene layers stacked together. 
Graphene layers placed together do not lie exactly one 
on top of each other but are shifted in such a way that 
only half of the carbon atoms have a neighbor in an- 
other layer and the other half are projected right into 
the middle of the graphene's "honeycomb cell". If the 
third layer aligns with the first (and the n + 2 layer with 
the n-th) then we arrive at the more stable arrangement 
of graphene layers known as Bcrnal (or AB) stacking, see 
Fig. [1^. However, this is not the only possible config- 
uration. One can imagine an alternative stacking when 
the third layer aligns neither the first nor the second but 
shifted with respect to both, see Fig. [TId. This arrange- 
ment is known as rhombohedral (or ABC) stacking and 
represents the main topic of this work. 

The very first studies of Bernal and rhombohedral 
graphitea^^— relying on the tight binding model have 
demonstrated the strong dependence of the band struc- 
ture on the stacking order. Later on the progress in nu- 
merical methods made it possible to refine the tight bind- 
ing model outcomes using ab initio calculations^^—. The 
seminal transport measurements on graphene^ have in- 
spired recent investigationsi^r— of the band structure in 
a few layer graphene having different stacking patterns. 
The band structure demonstrates a large variety of be- 
havior including a gapless spectrum, direct and indirect 
band gaps, and energetic overlap of the conduction and 
valence bands even though the number of layers has been 
limited by four In addition, the Landau levels and elec- 
tric field induced gap in a few layer graphene have been 
studied in RefsAS. ancU^ respectively. A very interest- 
ing effect of Landau-level crossing has been found very 
recently in trilayer grapheneJ^ It has also been shown 
that electron-electron interaction effects play an impor- 
tant role in a few layer graphene where charge carriers 
may exhibit a variety of broken symmetry statesJ^ 

The brief literature overview given above illustrates 
how rich the band structure of multilayer graphene 
(and the effects associated with) can be. Note, how- 




FIG. 1: (a) Bernal graphite represents graphene layers placed 
together in ...-AB-AB-AB-... stacking when the two lay- 
ers shifted with respect to each other by 0.142 nm alter- 
nate in bulk, (b) Rhombohedral graphite requires three non- 
equivalent graphene layers when the two layers are shifted 
with respect to the first one by 0.142 nm and 0.284 nm alter- 
nating in ...-ABC-ABC- ABC-... sequence. Dashed lines are 
shown for an eye guide. 

ever, that among all stacking possibilities, only the pure 
ABC arrangement maintains the sublattice pseudospin 
chiralityJi In the simplest case of negligible interlayer 
asymmetries and trigonal warp the simplified two-band 
ABC graphene model leads to the following Hamiltonian 
close to the neutrality point^i 

where vq w 10^ms~^ is the group velocity for carriers 
in single layer graphene, 71 is the hopping parameter, 
J/ = ± is the valley index, and N is the number of layers. 
As one can see from Eq. ([T]), the N- layer and (N-l-1)- 
layer graphene stacks differ, apart from the density of 
states, by only the winding number associated with the 
pseudospin orientation. The pseudospin texture in the 
momentum space associated with Hamiltonian ([T|) and 
shown in Fig. [2] is the main topic of this work. In what 
follows we focus on the influence of electron-electron ex- 
change interactions on the pseudospin texture and its 
detection by optical means. 

Note that the pseudospin lies in the xy-plane as long as 
its carrier remains in an eigenstate of Hq . The exchange 
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FIG. 2; Top view on the pseudospin texture for the conduc- 
tion band in chirally stacked multilayer graphene calculated 
from the non-interacting Hamiltonian ([l}. The pseudospin 
orientations in conduction and valence bands for a given mo- 
mentum are antiparallel. 



interactions can turn the pseudospin texture to the out- 
of-plane phase with the out-of-plane angle depending 
on the absolute value of the particle momentumi^iiSr^. 
This is due to the huge negative contribution to the 
Hartree-Fock ground state energy from the valence band 
(i. e. "antiparticle" states) which cannot be neglected in 
graphene because of the zero gap and conduction- valence 
band coupling via pseudospin. The broken symmetry 
states in multilayer graphene with chiral stacking have 
been recently studied by Zhang et aA^ using quite general 
arguments. Spontaneous symmetry breaking can occur 
in presence of exchange electron-electron interactions^*^. 
Here, we utilize a simplified model which includes ex- 
change electron-electron interactions but, at the same 
time, allows transparent half-analytical solution. Having 
this solution at hand we focus on the manifestation of 
such broken symmetry states in optical absorption mea- 
surements. 

Optical absorption via the direct interband optical 
transitions in monolayer graphene has been investigated 
irt2i and shown to be equal to the universal value ne'^ /he. 
In the presence of the electron-electron interactions the 
interband absorption can be substantially reduced or en- 
hanced as compared to its universal value ire^ /he just by 
switching the helicity of the circularly polarized light 
This effect is due to the peculiar pseudospin texture aris- 
ing from the interplay between pseudospin-momentum 
coupling and exchange interactions. To observe the pseu- 
dospin texture in chirally stacked multilayer graphene by 
optical means the photon energy must be much smaller 
than the bottom of the lowest split-off bands 71 = 0.4 eV. 
To give an example, a CH3OH 20 mW laser— with wave- 
length 118 /im (i.e. with the photon energy 10.5 meV) 
safely satisfies this condition. 



II. MODEL 

We start from the Coulomb exchange Hamiltonian for 
chiral carriers which is given by 



^oxch(l^) 



E 



1^ 



Ui^-^'l\x'k'.'){x'k'J (2) 



with ?7|k-k'| = 27re^/e|k — k'| and k' = ± being the 
band index with k = -I- for the conduction band. In 
order to consider the exchange Hamiltonian for any N 
on the equal footing we assume strictly two-dimensional 
Fourier transform for Coulomb potential. This is in 
contrast to^S where the interlayer distance d has been 
included in the screening multiplier. The intervalley 
overlap is assumed to be negligible, and the eigenstates 
-^cxch '^^'^ tie formulated as ^'^^(1") — 



of H"" = 



IxL) with spinors \xl+) = (cos f^, sin e'^^''^) 
= {niu^^-vcos^e^^'fY, and i&mp = ky/k^, 
where a non-zero out-of-plane pseudospin component 
corresponds to 1?^ ^ 7r/2. To diagonalize i?" the follow- 
ing j/-independent equation for dk must be satisfiedi^i^ 



(hvok) 



N 



■ COS -dk 



(-71)^ 

svadk' cosdk cos{Nip' ~ Nip)] — 
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Here the integration goes over the occupied states. Note 
that the conduction and valence states are entangled, 
and the latter cannot be disregarded even at positive 
Fermi energies assumed below. Thus, in order to eval- 
uate the integrals in Eq. ^ a momentum cut-off A is 
necessary. Most natural choice A = ji/hvQ corresponds 
to the energy scale 71 at which the split-off bands of bi- 
layer graphene become relevant and our two-band model 
no longer applies. Substituting x ~ k/A we arrive at 



N n 
X COSWfc 

a* 

2ir 1 

dip' 
kp/A 



. ,cosi?fc/ sini^fe — siniJfe/ cosiJfc COS A^iys' 
dx X — 



yj x'^ + X'^ — 2xx' COS p' 



(4) 



where a* = e'^/{ehvQ). The momentum cut-off is as- 
sumed to be much larger than the Fermi momentum kp, 
and, therefore, we can set the lower integral limit to zero. 
In this case our outcomes do not depend on the value of 
A. 

Besides a trivial solution with do — tt/2 independent of 
k, there are non-trivial ones t?i = 'd{k) and t?2 = Ti — d{k) 
with 'd{k) shown in Figs. [3^ and|4^ for different N and a* . 
The solutions 'da and ^1,2 represent to two phases with 
different total ground state energies E^^^ (-^'tot*) for the 
in-plane (out-of-plane) pseudospin phase. The difference 
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FIG. 3: Left panel: The pseudospin out-of-plane angle i9(fc) 
for substrate-free chirally stacked A^-layer graphene numeri- 
cally calculated from Eq. ((4]). Right panel: The total ground 
state energy difference ((5| between the in-plane and out-of- 
plane phases for different A'^. Increasing N makes the out- 
of-plane phase more preferable. The curve for = 1 
differs from the one given inJ^ since we have improved the 
precision of our calculations here. Note, that the Coulomb 
interactions are completely unscreened here. The screening 
is taken into account in Fig. |4] by introducing the effective 
dielectric constant e which relates to a* as a* = /{ehvo). 
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per volume is given by 
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XX 
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The ground state energy is the same for both valleys and 
spins, and, therefore, Eq. ([5]) contain gs — 2 and gy — 2 
for spin and valley degeneracy respectively. AiJtot has 
been evaluated numerically and the resulting AiJtot vs. 
N dependency is shovifn in Figs.[3)D and|3]3 for suspended 
and Si02-placed graphene respectively. One can see that 
strong electron-electron interactions with a* — 2 defi- 
nitely make the out-of-plane phase energetically prefer- 
able for > 1. For graphene placed on Si02 substrate 
the pseudospin out-of-plane phase is energetically prefer- 
able only for iV > 3. Note that the estimates of a* for 
clean monoatomic graphene flake vary from 2 (Ref<^) 
to 2.8 (Refi^) and, therefore, even monolayer graphene 
may get to the pseudospin out-of-plane phase. 

Note that it is possible to choose either the same or op- 
posite solutions for two valleys. The former choice breaks 
the parity invariance whereas the latter one does so with 
the time reversal symmetryjii^ In what follows we consider 
possible manifestations of these solutions in the optical 
absorption measurements on multilayer graphene. 



III. OPTICAL ABSORPTION 

We assume that an electromagnetic wave is propagat- 
ing along z axis perpendicular to graphene plane. It can 
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FIG. 4: Left panel: The pseudospin out-of-plane angle T}{k) 
for chirally stacked A^-layer graphene on Si02 numerically cal- 
culated from Eq. Q. Right panel: The total ground state 
energy difference ([SJ between the in-plane and out-of-plane 
phases for different A''. The difference is much larger in the 
case of substrate-free graphene, see Fig. |31 



be described by the electric field E = Eg exp{ikzZ — iujt), 
where kz is the wave vector, and uj is the radiation fre- 
quency. For our purpose it is more convenient to char- 
acterize the electromagnetic wave by its vector potential 
A = AoBxp^ikzZ — iiot), where Aq = icE/tu. The rela- 
tion between E and A is unambiguous since A is always 
in xy plane resulting in a gauge with divA = 0. The in- 
teraction Hamiltonian can be then written in the lowest 
order in A as 

jiyhvok)^ eN 
fc(-7i)^-i he 











(6) 



The total absorption P can be calculated as a ratio 
between the total electromagnetic power absorbed by 
graphene per unit square Wa — hujJ2v J j:^w{'k+,'k—) 
and the incident energy flux Wi — ui'^A^/Attc. The prob- 
ability to excite an electron from the valence band to an 
unoccupied state in the conduction band z/;(k-|-,k— ) can 
be calculated using Fermi's golden-rule 



Eh 



hiu). 



Because of normal incidence there is zero momentum 
transfer from photons to electrons. What is important 
is that the inter-band transition matrix elements of H^^^ 
turn out to be sensitive to the light polarization and pseu- 
dospin orientations in the initial and final states. In par- 
ticular, for linear polarization with = Acosiy9poi and 
Ay = A sin (y9poi the probability to excite an electron with 
a given momentum reads 

xS{Ek+ - Ek- - hw) 

X [l + cos^ i9fe - sin^ -dk cos{2Nip - 2Nippo\)] . (7) 

Note that w(k-|-,k— ) does not depend on the valley in- 
dex ly. The total optical absorption P is not sensitive to 
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the particular orientation of the polarization plane since 
the (/3poi-dependent term is integrated out in this case. 
In contrast, if we assume a circular polarization fulfilling 
Ax = ±iA/^/2, Ay — A/^/2 then the interband transi- 
tion probability for K-valley can be written as 



/, , M 47r f eN / 

w;(k+,k-)|,=+ = — i—A 



2Ni,2N-2 



(-71) 



2N-2 



^ ^ cos^ 4 ^^^'■+ " " 



(8) 



Here, the multipliers sin'*(i?i:/2) and cos'*(i9fe/2) are for 
two opposite helicities of light, and for K'-valley they are 
interchanged. If the out-of-plane pseudospin polarization 
is chosen to be opposite in two valleys then this helic- 
ity dependence survives the integration over momentum 
and summation over valley index resulting in the helicity- 
sensitive total absorption which reads 



2N 



P 



A 



16N^ Tre^ {hv„) 

~luL (-71)2^-2 
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X I dkk 





2N-1 



sm 



cos 



The absorption strongly depends on helicity as long as 
the radiation frequency is much smaller than the band 
split-off energy. This regime corresponds to the exci- 
tation of electrons with comparatively small momenta, 
where the angle -dk is close to zero, see Fig. |3Jd and 
Fig. ib. 

The experiment proposed here is similar to the one 
discussed recently in Ref.^^ suggesting to detect broken 
symmetry states in graphene placed on substrate by polar 
Kerr rotation measurements, i.e. by analyzing the ellip- 
tical polarization of reflected radiation. In the present 
paper, however, we propose to reach the same goal by 
comparing the opacity for two orthogonal light polariza- 
tions, which appears to be the simpler strategy. 

Alternatively, one can characterize broken-symmetry 
states in terms of the topological Chern numbers^! given 
by 



In J 



(10) 



with A!^ = i{Xk±\^k\Xk±)- This approach is routinely 
used in the theory of topological insulators^^ and orig- 
inated from the description of quantized Hall conduc- 
tances. 

A non-trivial band topology has also been found in 
graphene^^^ where the band gap was opened via spin- 
orbit coupling characterized just by a constant (i.e. wave 
vector independent) mass term rather than by exchange 
interactions described by a more complicated Hamilto- 
nian 

Note, that in our case here x Aj. contains a term 
proportional to S(k) which occurs due to the singularity 



in A'j. (X but does not contribute to as long as 
i?fc=Oatfc— >0. Thus, for the conduction band electrons 
we have 



Mh sin 'dk 



(11) 



Assuming the out-of-plane pseudospin polarization to be 
opposite in two valleys it follows that the total Chern 
number for conduction electrons just equals N . On the 
other hand, the total Chern number is zero as long as the 
out-of-plane pseudospin pseudospin component is same 
in both valleys. Thus, the two out-of-plane solutions 
are topologically different even though they correspond 
to the band gap of the same size. The topologically non- 
trivial case with n 7^ corresponds to broken time rever- 
sal symmetry leading to the existence of a zero-field Hall 
currenti^. 

The only question is whether the time reversal bro- 
ken states really occur in clean graphene samples. 
More in depth analysis performed by Nandkishore and 
Levito v^^i'^^ suggests that this is the case at least for bi- 
layer graphene. On the other hand Jung et alM^ demon- 
strate that the intervalley exchange coupling favors par- 
allel pseudospin polarization in opposite valleys breaking 
the parity invariancc. In this case the total absorption 
does not depend on the radiation helicity but the two val- 
leys are occupied differently by the photoexcited carriers 
which might be useful effect for valleytronics^. 

In the in-plane phase with d = i: /2 the total absorption 
does not depend on light polarization, and in the non- 
interacting limit it equals to 



P ^ N 



ne 

he 



(12) 



At N -- 
expect ed^i. 



1 it acquires the universal value ^ 



he 



IV. CONCLUSION 

We have demonstrated that the polarization-sensitive 
optical absorption predicted iriF^ for single layer graphene 
can be found in chirally stacked carbon multilayers 
within a broader range of parameters. This is due 
to the enhancement of the interaction effects in multi- 
layer graphene with larger number of layers. The con- 
ditions necessary for the observation of the polarization- 
dependent optical absorption can be summarized as fol- 
lows, (i) Graphene samples must be prepared as cleaner 
as possible. It would be better to utilize suspended sam- 
ples in order to isolate the interacting electrons from the 
environment with large relative permittivity, (ii) The ex- 
change interactions must break the time reversal invari- 
ancc rather than the parity. This situation corresponds 
to the topologically non-trivial state with the Chern num- 
ber equal to N. (iii) Finally, the low frequency radiation 
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is necessary to excite the electrons with smaUcr momenta 
having larger out-of-plane pseudospin components and to 
preclude the influence of split-off bands. 

In addition, we would like to mention the possible ex- 
istence of a zero-field Hall current in the slightly doped 
graphene samples where the time reversal invariance is 
broken as described above. 
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